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NOTATION 


inertia  matrix 

^mn 

elements  of  the  inertia  matrix 

b 

damping  coefficient  for  the  hinges 

^mn 

elements  of  the  stiffness  matrix 

c 

spring  stiffness  in  hinges 

denominator  in  reduced  frequency  expression 

DR 

divergence  region 

FB 

flutter  boundary 

G 

m ' g Jl/c 

g 

acceleration  of  gravity 

k 1 , k 2 

stiffness  of  horizontal  springs 

£ 

length  of  a bar  in  the  double  pendulum 

m ^ 2 

concentrated  masses 

m * 

a mass  parameter 

m , n , N 

positive  integer 

N] 

numerator  in  reduced  frequency  expression 

P 

magnitude  of  the  applied  circulatory  force 

Pj 

coefficients  in  equation  (29),  j=l,2,3,4 

0 

^d’^d 

Pl/c 

critical  flutter  loads  in  the  presence  of  arbitrary  and 
very  slight  damping,  respectively 

Qg  »Qg 
1 2 

« 

AQ 

critical  loads  of  divergence  of  the  first  and  second 
types,  respectively 

difference  of  critical  loads 

% 

generalized  coordinates 

SR 

stable  region 

4- 

U 

time 

iii 
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± 


constants  of  integration 
a tangency  coefficient 

Aj  characteristic  determinants,  j=1,2,3 

n b/ac,  dimensionless  damping  parameter 

K j £^/c , i = 1 »2 

X characteristic  exponent 

X|^,Xj  real  and  imaginary  parts,  respectively,  of  X 
^2  a dimensionless  mass  parameter 

a (m'£^/c)^^^ 

T dimensionless  time 

(})■],(}) 2 angles  specifying  the  configuration  of  the  double 

pendulum 


IV 


INTRODUCTION 


The  differential  equations  of  motion  of  a non-di ssi pati ve , 
non-gyroscopi c system  whose  configuration  is  specified  by  N 
generalized  coordinates  R]  >^2  ’ * * * generally  of  the  form 

^ (^mn^n  ^mn^n^  “ m , n = l , 2 , • • • , N , (1) 

n = l 

"2  2 

where  q^  = d q^^/dt  and  and  denote  respectively  the  elements 
of  the  inertia  and  stiffness  matrices  for  the  system.  Upon  per- 
formance of  equivalent  operations  such  as  multiplying  the  individual 
equations  of  motion  by  constants  and  adding  the  resulting  equations 
to  one  another,  the  form  of  the  set  of  equations  in  (1)  can  be 
modified.  Ku^  calls  the  system  in  (1)  a purely  dynamic  system  if 
the  number  of  second  order  differential  equations  in  (1)  remains 
unchanged  under  the  application  of  any  sequence  of  equivalent 
operations.  If,  however,  apposite  equivalent  operations  degenerate 
one  or  more  of  the  second  order  differential  equations  into  purely 
algebraic  equations,  Ku  calls  such  a system  a qua s i - dy nami c system. 
Indeed,  if  the  matrix  A=(Aj^p)  is  non-singular,  the  system  is  purely 
dynamic,  whereas  if  det  (A)  = 0,  then  the  system  is  quasi -dynamic . 

In  subsequent  investigations,  Ku^’^’^  determined  that  a quasi- 
dynamic system  consisting  of  a double  pendulum  composed  of  two 
elastically  hinged  weightless,  rigid  bars  carrying  only  a single 
concentrated  mass  and  subjected  to  a tangential  force  applied  at 
its  free  end  could  lose  stability  only  by  divergence.  Furthermore, 

*In  the  present  investigation  the  term  "divergence"  means  that 
characteristic  displacements  of  the  system  grow  monotonically  and 
exponentially  with  time. 
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he  showed  that  the  onset  of  divergence  for  this  particular  quasi- 
dynamic system  Is  characterized  by  Its  sole  natural  frequency  of 
vibration  becoming  Infinite.  But,  as  will  be  demonstrated  herein, 
these  observations  shed  light  on  only  one  aspect  of  the  possible 
stability  behavior  that ' quasi -dynamic  systems  are  capable  of 
exhibiting.  In  the  present  study,  Ku ' s applied  tangential  force 
Is  replaced  by  a force  whose  orientation  Is  specified  by  a cer- 
tain tangency  coefficient.  It  will  be  demonstrated  here  that 
such  a system  can  become  unstable  only  through  divergence  with 
the  onset  of  Instability  being  signaled  either  by  the  vanishing 
of  the  frequency  or  by  the  frequency  becoming  Infinite.  In 
addition,  it  will  be  shown  that  stability  maps  for  a quas 1 -dynami c 
system  often  possess  multiple  regions  of  stability  and  divergence. 

It  Is  also  of  interest  to  investigate  the  influence  of  damping 
on  the  stability  of  discrete  systems  that  have  singular  Inertia 
matrices.  For  lack  of  more  suitable  terminology,  the  term  "quasi- 
dynamic", which  was  applied  originally  only  to  non-di ssipatl ve 
systems,  will  be  retained  throughout  the  present  discussion,  where 
it  will  be  understood  that  any  non-gyroscopic  system  for  which 
det(A)=0  may  be  called  quasi -dynamic  independent  of  the  presence 

a. 

or  absence  of  damping  that  Is  linearly  proportional  to  velocity. 

The  results  reported  here  for  the  non-conservati vely  loaded  double 
pendulum  having  a singular  Inertia  matrix  and  viscoelastically 
hinged  joints  indicate,  as  one  should  expect,  that  Internal  damping 
has  no  effect  whatsoever  on  the  value  of  the  critical  load  of 
divergence,  whose  onset  is  characterized  by  a vanishing  frequency. 
In  marked  contrast,  however,  it  will  be  shown  that,  for  certain 
values  of  the  tangency  coefficient,  the  damped  system  loses 
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stability  by  flutter,  whereas  its  non-di s s i pa t i ve  counterpart 
becomes  unstable  by  a form  of  divergence  whose  onset  is 
characterized  by  a frequency  becoming  infinite.  Furthermore, 
for  arbitrarily  small,  positive  internal  damping,  the  value 
of  the  critical  flutter  load  is  observed  to  be  less  than  the 
value  of  the  critical  divergence  load  for  the  corresponding 
undamped  quasi -dynamic  system. 

EQUATIONS  OF  MOTION 

The  double  pendulum  depicted  in  Figure  1 consists  of  two 
rigid,  weightless  bars  of  length  A concentrated  mass  m2 

is  attached  to  the  upper  bar  at  the  point  B,  and  the  weight 
m2g  of  this  mass,  where  g denotes  the  acceleration  of  gravity, 
acts  vertically  downward.  Frequently,  in  studies  of  the 
stability  of  a double  pendulum  subjected  to  a non-conservative 
force,  a second  mass  m-|  is  affixed  to  the  lower  bar,  often  at 
the  point  B.  In  this  study,  however,  it  will  be  assumed  that 
m.|=0,  which  implies  immediately  that  the  inertia  matrix  of  the 
system  will  be  singular.  Thus,  the  system  will  herein  be 
called  quas i -dynami c . At  the  joints  0^  and  A , the  double 
pendulum  has  viscoelastic  hinges  that  exert  linear  restoring 
moments  c<j>i+b^i  and  c((})2-<}>i)  + b{^2“^’T)>  where  c,  b>0  are  a 
stiffness  coefficient  and  an  internal  damping  coefficient, 
respectively.  The  system's  configuration  is  specified  by  the 
(assumed  small)  angles  4>i  and  4>2  formed  between  the  vertical 
and  each  of  the  bars.  The  points  A and  B are  supported  by 
horizontal  linear  springs  of  stiffness  k-]  and  k^,  while  a 
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compressive  force  P oriented  at  an  angle  a(fi2  relative  to  the 
vertical,  where  a denotes  the  tangency  coefficient,  is  applied 
at  the  end  B. 

The  equations  of  motion  for  this  system  may  be  shown  to  be 

+ 2b  + (2c-PJl+k^  Jl^+k2il^-m2gil)4)i  + - 

- b ^2  ■ (c-k2A^-aPJl)4)2  = 0,  (2) 

- b - (c-k2i.^)4>i  + + b ^2  + tc  - 

- PJl(l-a)  + k2il^-m2gJl]4>2  = 0*  (3) 

where  ^^=d4>^/dt,  etc.  To  facilitate  the  analysis  of  these 
equations,  it  is  expedient  to  render  them  in  dimensionless 
form.  Defining  m2  = y2’’’' » Kj  = kji-3/c^  6=m'gA/c,  n=b/ac, 

and  t = OT,  where  j = l,2  and  a^=m'il^/c,  one  finds  that  (2)  and 
(3)  become 

^ ••  • 

P24>i  + 2n  ^1  + (2-Q  + k-| +K2-y2^^1  ^2*^2  " ‘^’2  “ 

- (l-aQ-K2)4>2  “ 

••  • "*  • 

y2‘J’i  " n <j)-|  - CI-K2H1  + y24>2  + n 4>2  + 

+ [1+Q(a-1)  + <2  “ ^2^^'**2  ^ 

where  now  (j>i=d4)i/dT  etc.  The  elements  of  the  inertia  matrix  A 
associated  with  (4)  and  (5)  are,  clearly,  A^^=y2  for  m,n=l,2. 
Hence,  ^ is  singular. 
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THE  NON-DISSIPATIVE  SYSTEM 


If  no  damping  is  present  in  the  hinges,  then  ni“n2~0  and 
(4)  and  (5)  reduce  to 

y2<l>T''(2-Q+Ki ^2^2  " n-aQ-K2^*f’2  “ 


y2‘^T  “'^2  '*'^2^2'''^^ 

If  (7)  is  subtracted  from  (6),  one  finds  that  the  inertia  terms 
Vi2^]  and  ^2*^2  cancel  and  that  a purely  algebraic  relationship, 
namely , 


4>2  = (3  + Ki -y2G-Q)<t»i/(2-y2^"Q) » (8) 

exists  between  4>i  and  <p2-  Ku^  calls  a relationship  of  the  type 
in  (8)  an  i nternal  constraint. 

Inserting  (8)  into  either  of  the  differential  equations  in 
(6)  or  (7),  one  obtains  the  following  reduced  equation  of  motion 
for  the  quas i -dynami cal  system: 

+ 0)^%^  =0,  (9) 

where 

= N^/D^  (10) 

with 

N-|  = (l-a)Q^  - [ ( 1 -a)  ( 3 + <i -U2G)  - P2^  "*■  "*■ 

+ 1 + <1  + <2 ( ^■*■'^1  ) "Vi2^  ^ ^^*^2  ^ ^ {vi2^^^ 

(11) 
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and 

= y2(5  + i<i-2U2G-2Q).  (12) 

It  is  obvious  from  (9)  that  the  motion  of  the  system  is  characterized 
by  a single  natural  frequency  to-i  . Consequently,  the  occurrence  of 
flutter,  which  ensues  upon  the  coalescence  of  two  natural  frequencies 
when  the  loading  conditions  are  suitable,  is  impossible  in  the  quasi- 
dynamic system  considered  here  even  if  its  purely  dynamic  counterpart 
is  flutter  prone  for  a given  circulatory  force.  Therefore,  the  state 
of  stability  of  the  given  quasi-dynamic  system  can  now  be  determined 
from  the  algebraic  sign  of 

If  u^>0,  it  is  evident  from  (9)  that  the  motion  of  the  system  is 

stable  since  it  consists  of  a bounded  simple  harmonic  oscillation. 

2 

On  the  other  hand,  if  toi<0,  the  motion  grows  exponentially  in  time; 
hence  the  system  is  unstable  through  divergence.  In  view  of  (10), 
it  can  now  be  stated  the  given  quasi -dynamic  system  will  be  stable 
whenever  N-|  and  D-|  have  the  same  algebraic  sign  and  unstable  by 
divergence  whenever  N-|  and  D-|  have  opposite  signs.  The  boundaries 
of  the  regions  of  stability  and  instability  in  stability  map  can  be 
obtained  from  the  conditions 

= 0.  (13) 

For  certain  combinations  of  parameters,  however,  it  happens  that  N-| 
cannot  vanish.  In  such  cases,  the  stability  boundary  is  determined 
from  only  one  of  the  conditions  in  (13),  namely, 

= 0.  (14) 
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If  O'!  is  persistent  in  sign  while  N-|  changes  sign  as  Q is 
increased,  the  system  becomes  unstable  by  divergence  and  the 
critical  divergence  load  may  be  computed  from  the  condition 
of  vanishing  frequency,  i.e.,  o)^  = 0,  which  implies  that  N^=0. 
Instability  whose  onset  is  characterized  by  the  vanishing  of 
will  be  called  here  divergence  of  the  first  type . It  may 
also  be  remarked  that  a rather  unusual  feature  of  the  present 

quasi -dynamic  system  is  that  its  motion  can  become  divergent 

2 7 

even  though  does  not  vanish  when  the  sign  of  is  changed 

from  positive  to  negative  as  the  value  of  the  load  parameter 

Q is  increased  up  to  and  beyond  its  critical  value.  Specifically 

suppose  that,  for  a given  value  of  Q,  both  and  D-|  are  positive 

If,  upon  increasing  Q,  N.j  remains  positive  while  D.j  becomes 

negative,  then  N-j  and  D.j  have  opposite  signs.  Thus,  the 

motion  of  the  system  must  be  divergent.  Clearly,  the  value 

of  does  not  vanish  at  the  onset  of  instability,  but  it  does 

become  infinite.  This  phenomenon  is  termed  divergence  of  the 

second  type . If  both  N.j  and  D-|  vanish  for  some  value  of  Q, 

then  the  value  of  must  be  obtained  with  the  aid  of 

L'Hospital's  rule.  In  references  [2]-[4],  Ku  discussed  only 

divergence  of  the  second  type  since  the  system  that  he  studied 

was  subjected  to  a tangential  force  which  cannot  produce 

divergence  of  the  first  type. 
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SPECIAL  CASES  OF  THE  NON-D I SS I PAT  I V E SYSTEM 

The  following  special  systems  will  be  considered  in  some  detail: 
(a)  G=k^=K2=0,  the  effect  of  weight  is  Ignored  and  the  two  horizontal 
springs  are  absent;  (b)  k^=<2=0,  the  effect  of  weight  is  retained  but 
the  horizontal  springs  are  absent;  (c)  G = k-|=0,  the  spring  at  B in 
Figure  1 is  included  in  the  system  but  the  effect  of  weight  is  ignored 
and  the  spring  at  A is  absent;  (d)  G=K2=0,  the  spring  at  A is  present 
in  the  system,  whereas  the  effect  of  weight  is  ignored  and  the  spring 
at  B is  absent. 

For  Case  a,  (11)  and  (12)  become 

= (l-a)Q2  - 3(l-a)Q  + 1,  (15) 


= P2(5-2Q), 

so  that  (10)  can  be  expressed  as 


(16) 


= [(l-a)Q^  - 3(l-a)Q+l]/p2(5-2Q).  (17) 


Let  Qg^  and  Qg^  denote  the  critical  loads  of  divergence  of  the 
first  and  second  types,  respectively.  Then,  to  determine  the 
stability  boundaries,  the  conditions  in  (13)  in  conjunction  with 
(15)  and  (16)  yield 

(l-a)Q^^  - 3(l-a)Q„  +1=0,  (18) 

®1  ®1 

and 


Q„  = 5/2.  (19) 

®2 
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Equation  (18)  was  first  studied  by  Herrmann  and  Bungay^,  and  (19) 
was  obtained  in  references  [2]-[4]  in  the  special  case  of  a=l . 

To  identify  the  regions  of  stability  and  instability  for  the 

given  system,  one  need  only  examine  the  eigencurves  obtained  from 

(17).  Representative  eigencurves  prepared  with  vi2“^  (with  no  loss 

of  generality  are  shown  in  Figure  2 for  four  values  of  the  tangency 

coefficient  a.  These  eigencurves  consist  of  two  branches  which 

asymptotically  approach  the  horizontal  line  Q=5/2  as  For  a 

given  value  of  a,  these  eigencurves  reveal  the  possibility  of 

multiple  intervals  of  stability  and  instability  relative  to  the  load 

parameter  Q.  The  case  of  a=l/5  is  exceptional,  however,  because  (17) 

assumes  the  rather  simple  form  = — (1-2Q).  Hence,  it  is  clear  that 
2 

0)1  decreases  linearly  in  Q and  the  system  loses  stability  by  divergence 
of  the  first  type  at  Q-1/2.  A stability  map  in  the  Qa-plane  that  is 
obtained  from  (18)  and  (19)  appears  in  Figure  3.  The  various  stable 
regions  (SR)  and  divergence  regions  (DR)  have  been  identified  with 
the  aid  of  the  eigencurves  shown  in  Figure  2. 

As  a second  example,  consider  the  double  pendulum  in  a gravity 
field  (g40).  case  b.  In  the  discussion  that  follows,  both  positive 
and  negative  values  of  G,  referring  to  either  a vertically  upright 
and  a vertically  hanging  double  pendulum,  respectively,  will  be 
considered.  Under  the  assumptions  that  <i  = <2  = 0 and  ^2  " 1 ♦ 
(11)-(13)  yield 

(l-a)QeJ  - [(l-a)(3-G)-G]Q^  + 1 - 3G  + G^  = 0,  (20) 

Q = \ (5-2G),  (21) 

2 2 
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the  latter  being  a straight  line  having  a slope  of  -1  In  the 
QG-plane.  When  a=l , (20)  leads  to 

Q-  = - (1-3G  + g2)/G.  (22) 

In  Figure  4,  stability  maps  In  the  QG-plane,  as  obtained 
from  (20)-(22)  are  presented  for  a=0,  3/4,  and  1.  In  Figure 
4a,  all  the  stability  boundaries  are  parallel  straight  lines, 
namely,  that  given  In  (21)  as  well  as  Q = ^ (3  +/5’  - 2G), 
which  are  obtained  from  (20)  when  a=0.  Upon  examining  the 
appropriate  elgencurves  for  the  given  system,  one  can  conclude 
that  the  stability  maps  shown  In  Figure  4 possess  multiple 
regions  of  stability  and  divergence,  whose  boundaries  depend 
strongly  upon  the  weight  parameter  G and  the  tangency 
coefficient  a. 

For  a non-conservatl vel y loaded  double  pendulum  whose 
motion  (case  c)  is  constrained  by  the  presence  of  a horizontal 
linear  restoring  spring  acting  at  its  upper  end  B,  as  depicted 
In  Figure  1,  the  stability  boundaries  are  derivable  from 

(l-a)Qe^  - [3(1-a)  + 2k  ]Q  + 1 + 5k.  = 0,  (23) 

Q.  = 5/2.  (24) 

2 

Equation  (23)  was  previously  given  by  Sundararajan^,  and  (24)  is 
Identical  to  (19)  which  was  obtained  for  the  system  in  which  the 
end  support  was  absent.  The  stability  boundaries  found  from  (23)- 
(24)  are  plotted  In  Figure  5 for  a=0,  1/4,  1,  and  3/2.  It  is 
evident  that  the  forms  of  these  curves  are  very  significantly 
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Influenced  by  the  values  of  the  elastic  constraint  parameter  <2 
and  the  tangency  coefficient  a. 

In  the  case  just  described,  the  stability  boundary  obtained 
from  (24)  Is  merely  a horizontal  straight  line  which  Is  Independent 
of  <2*  contrast.  If  the  spring  Is  removed  from  the  end  B (see 
Figure  1)  and  Is  Instead  attached  horizontally  at  the  hinge  A In 
the  double  pendulum,  (11)-(13)  yield,  for  case  d, 

(l-a)Q  ^ - (l-a)(3  + <jQ_  + 1 + <,  = 0 , (25) 

®1  1 ®1  ' 

and 

Qe^  = \ (5+<^).  (26) 

It  should  be  noted  that  when  a=l  , (10)  becomes 
= (l+<^)/y2(5  + <i-2Q). 

Consequently,  In  the  case  of  a tangential  force,  stability  Is 
lost  only  through  divergence  of  the  second  type  with  the 
critical  divergence  load  being  obtainable  from  (24). 

The  stability  boundaries  obtained  from  (25)  and  (26)  have 
been  plotted  in  the  stability  maps  In  Figure  6 for  five  values 
of  a,  where  multiple  regions  of  stability  and  divergence,  whose 
extents  depend  strongly  upon  a,  are  shown.  The  secondary  (upper) 
stability  regions  appearing  In  Figure  6 change  their  character  in 
the  sense  of  separating  from  each  other  as  the  value  of  a Is 
Increased  through  and  beyond  a=l/2.  Indeed,  for  a=l/2,  (25) 
assumes  the  form 

Qg^  - (3+<i )Qg^  + 2(l+<i ) = 0, 
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whence 


Qg  ■ ^ (3+ic,i|K,-l  I).  (27) 

Selecting  the  negative  sign  in  (27),  one  finds 

if  0<K^<1, 

i f , 

whereas,  upon  choosing  the  positive  sign, 

2 if  0<<^<1, 

1 +<i  i f 1 _<K-|  . 

Therefore,  these  two  expressions  lead,  in  essence,  to  a pair  of 
intersecting  straight  lines,  as  shown  in  Figure  6c,  which  form 
the  boundaries  of  the  secondary  region  of  stability. 

THE  DISSIPATIVE  SYSTEM 

A solution  of  the  set  of  homogeneous  equations  (4)  and  (5), 
in  which  the  damping  terms  have  been  retained,  will  be  sought 
in  the  form 

<j>^(x)  = X^e^'^,  n=l,2,  (28) 

where  the  X^'s  are  constants  and  the  characteristic  exponent  X is, 

1/2 

in  general,  a complex  number,  i.e.,  X=Xj^+iXj,  where  i = (-  l)  ' . 


1 +Ki 


12 


Substitution  of  (28)  into  (4)  and  (5)  leads  to  a set  of  linear 
homogeneous  equations  from  which  is  derived  in  the  familiar 
manner  the  frequency  equation 


(29) 


where 


= n[2  + <i +5<2-3u2G"’'3Q  (a-1  ) ] t 


(30) 


P4  “ 1 + <1  + 5<2  + <1  <2  ■ ^2®  ( ■*'^'^2  ^ ^ 

+ Q[(3  + <^  )(cx-l  )-2<2  + (2-a)u2G]  + Q^(l-a). 


In  view  of  the  form  of  equation  (28),  it  is  clear  that  the 

motion  of  the  system  will  be  stable  if  all  three  values  of  X 

obtained  from  (29)  are  such  that  X[^<0.  On  the  other  hand,  if 

X >0  and  X.=0  for  at  least  one  of  the  three  characteristic 
R 1 

roots,  the  system  will  be  unstable  by  divergence,  whereas, 

if  ^[^>0  and  Xj^O,  for  at  least  one  of  the  roots  of  (29),  the 

motion  consists  of  an  oscillation  with  an  exponentially 

increasing  amplitude,  and  the  equilibrium  configuration  of 

the  system  is  said  to  be  unstable  by  flutter. 

The  Routh-Hurwi tz  criterion  states  that  the  necessary  and 

sufficient  conditions  for  all  the  roots  X.,  j=l,2,3,  of  (29) 

J 

to  have  negative  real  parts  is  that  p-|>0  and  all  the  determinants 
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^1  “ P2 » 


A 


2 


P2  Pi 

P4  P3 


P2  Pi  0 

P4  P3  P2 

0 0 0 


be  positive , i . e . , 

p^>0,  P2>0,  P2P3-PiP4>0,  P4>0. 


(31) 


since  A^=P4A2.  In  view  of  (30),  the  inequalities  in  (31) 
become 


p^  = 5y2n>0, 

Pg  = ^2^^  *^1  ■ ^^2^  ■ 

P2P3-P1P4  = n{y2t(2-y2*^^^  + (l+<l)^]  + n^[2  + Ki+5iC2-3y2G]- 
- [4  + 2k^ -y^G- (1 -a)  (2k^+P2G-3ti^)]Q  + 

+ y2(l-ct)Q^>>0.  (32) 

P4  = 1 + <1  + 5<2  + <1  <2  ■ y2^(^'*’'^l''’^'^2^''’(y2^^^  " 

- [(3+K^)(l-a)  + 2k2  - (2-a)y2G]Q  + (l-a)Q^>0. 

In  the  present  case,  the  critical  flutter  load  Q^j  is  determined 
from  the  condition  P2P3"PiP4=0.  which,  in  view  of  (32),  becomes  as 
long  as  n>0 

(l-a)Qj  - [4+2<i -G- (1 -a) (2Ki+G-3n^)]Qj  + (2-G)^  + 

+ (1+Kl)^  + ti^[2ki +5k2“3G]  = 0,  (33) 

upon  setting  P2  = l . with  no  loss  of  generality  for  present  purposes. 
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In  the  case  of  a tangential  force  (a=l),  (33)  yields 

Qd  = C(2-G)2+(1+<^)2+^2(2+^^+5^^.3G)]/(4+2k^-G).  (34) 

SPECIAL  CASES  OF  THE  DISSIPATIVE  SYSTEM 

In  analogy  with  the  special  cases  (a)-(d)  of  the  non-di ss i pa ti ve 
system  discussed  in  an  earlier  section,  attention  may  now  be  turned 
to  the  corresponding  situations  for  the  dissipative  system,  namely. 

Case  a.  G = <1  = <2  “ 0- 


(l-a)Qd  ■ [4+3n^(l -a)]Qj+5+2n^  = 0; 


(35) 


Case  b.  " *^2  " 


(l-a)Qj  - [4-G-(l-a)(G-3Ti^)]Qjj  + 1 + 
+ (2-G)^  + n^(2-3G)  = 0; 


(36) 


Case  c.  K,  = G = 0: 


(l-a)Qj  - [4+3n2(l-a)]Q^  + 5 + h^(2+5k2)  = 0;  (37) 


Case  d.  <2  " ^ " 3- 


(l-a)Qj  - [4  + 2a<i  + 3n^(l-a)]Qj  + 4 ' + 
+ (l+<,  )^  + ri^(2+K.  ) = 0. 


(38) 
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In  the  event  that  the  applied  load  is  a tangential  force  (a=l), 


(35)-(38)  yield 


= (5  + 2n^)/4. 


(39) 


= [1  + (2-G)2  + n^(2-3G)]/(4-G). 


(40) 


Q(1  = [5  + 0^(2  + 5<2)]/4. 


(41) 


Q(1  = [4  + (1  + <i)^  + n^(2  + <t)]/2(2  + <i).  (42) 


respectively. 

In  Figure  7,  the  stability  map,  i.e.,  a plot  of  the  variation 
of  the  critical  load  Q versus  the  tangency  coefficient  a,  has  been 
prepared  for  the  case  of  a double  pendulum  for  which  G = k-j  = <2  “ 
case  a.  The  solid  curves  form  the  divergence  boundaries  (OB) 
associated  with  divergence  of  the  first  type,  and  these  boundaries 
are  completely  independent  of  the  value  of  the  damping  parameter  n. 
The  dot-dash  curve,  associated  with  instability  due  to  divergence 
of  the  second  type,  also  represents  a divergence  boundary  which, 
however,  is  obtained  only  in  the  complete  absence  of  damping  (n=0). 
If  damping  is  present  (n>0)  in  the  system,  divergence  of  the  second 
type  no  longer  occurs,  but  now  instability  due  to  flutter  appears 
in  its  place.  The  dash-dash  curve,  which  provides  the  flutter 
boundary  (FB),  is  obtained  from 


(l-a)Q*^  - 4Q*  + 5 = 0 


(a+1) 


(43) 
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which  results  from  (43)  in  the  case  of  slight  damping,  i.e., 

0<n<<l . or  from  Q^=5/4  when  a=l  which  is  found  from  (39)  upon 
neglecting  n- 

If  n=0,  the  stability  region  (SR)  is  bounded  by  the  solid 
and  dot-dash  curves.  For  the  slightly  damped  system  (0<n<<l). 
the  solid  curves  are  still  pertinent,  but  now  the  dot-dash 
curve,  corresponding  to  divergence  of  the  second  type,  is  not 
since  the  dot-dash  divergence  boundary  must  be  replaced  by 
the  dash-dash  flutter  boundary.  Upon  examination  of  (18)  and 
Figure  7,  it  is  evident  that  flutter  is  possible  only  when 
a>5/9  and  n>0.  If  a<5/9,  the  system  becomes  unstable  by 
divergence  of  the  first  type.  Furthermore,  it  is  obvious 
that  Q(j<Qe  • 

In  Case  b,  (40)  becomes 

(l-a)Qj^  - (4-G)Q*  + 1 + (2-G)2  = 0 

when  0<n<<! . Stability  maps  in  the  QG-plane  are  shown  in  Figure  8 
for  a=3/4  and  a=l . In  both  cases,  the  stability  boundary  consists 
of  a flutter  boundary  and  a divergence  boundary  which  intersect  at 
a point  in  the  first  quadrant. 

If  the  upper  end  of  the  spring  is  supported  by  a spring  (k2"^0)» 
then,  with  k.|=G=0,  (41)  is  appropriate,  and  for  very  slight  damping 
this  equation  becomes  identical  to  (43).  Therefore,  in  the  presence 
of  slight  internal  damping,  the  value  of  the  critical  flutter  load 
is  independent  of  the  spring  stiffness  parameter  k^.-  The  flutter 
boundaries,  consequently,  are  identified  by  the  horizontal  dash-dash 
lines  in  Figures  9a-c,  which  have  been  prepared  for  a=l/4,  1,  and  3/2, 
respectively . 
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Finally,  the  flutter  boundary  for  a slightly  damped  double 
pendulum  supported  by  a spring  at  the  central  hinge  may  be 
determined  from 

(l-cx)Qd^  - 2(2+aK^)Q*  + 4 + (1+k^)^  = 0 

which  is  obtained  from  (42)  upon  neglecting  the  terms  containing 
n-  Two  stability  maps  in  the  Q<^-plane  are  presented  in  Figure  10 
for  a=0.504  and  a=3/5. 


SOME  COMPARISONS 

In  the  case  of  a tangential  force,  a=l , (34)  becomes 

Q*  = [(2-G)2  + (1+<^)2]/(4+2k^-G)  (44) 

for  very  slight  damping,  i.e.,  0<n<<! • For  the  corresponding 
undamped  system,  one  has 

Qe^  = ^ (5  + <1  - 2G),  (45) 

upon  setting  D-|=0  and  P2~^  (12).  The  difference  of  these 

quantities  is  easily  shown  to  be 


AQ  = Qe2  - Qd 

= 5(1+k^  )(2-G)/2(4  + 2k.|-G). 


(46) 


If  G=0,  then  (46)  becomes  AQ  = 5(1+k.|  )/2(2+k^  ) , which  is  always 
positive  for  k-|^( 
then  (46)  yields 


positive  for  k-|^0.  Thus,  Qd^^Qg  • other  hand,  k-|=0. 


AQ  = 5(2-G)/2(4-G), 


(47) 
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★ 

which  is  meaningful  only  for  Gf2/3.  Thus,  it  follows  from  (47) 

that  AQ^l  for  G£2/3.  Consequently,  Qj<Qp  . Therefore,  one  may 

d 

conclude,  in  the  case  of  a tangential  force,  that  the  value  of 
the  critical  flutter  load  for  the  slightly  damped  quasi -dynami c 
system  is  always  less  than  the  value  of  the  critical  load  of 
divergence  of  the  second  type  for  its  undamped  counterpart. 

In  the  other  extreme  of  very  large  damping,  it  follows  from 
(34)  for  a tangential  force  that  an  approximation  for  the 
critical  flutter  load  is 

Qjj  n^(2  + K-|+5K2-3G)/(4  + 2K.| -G)  , (48) 

from  which  it  is  clear  that  the  value  of  is  essentially 
proportional  to  the  square  of  the  internal  damping  parameter 
n.  In  contrast,  one  can  show  that  the  two  solutions  of  (33) 

p 

behave  as  Q^'^3n  and 

Qd  (2  + K^+5<2-3G)/3(l-a)  (49) 

as  provided  that  a=fl  . Equation  (49)  yields  the  critical 

flutter  load  for  very  large  n,  and,  in  addition,  it  also  yields 
the  maximum  critical  value  of  when  Qj  is  treated  as  a function 
of  q. 


*This  value  of  G is  obtained  from  3G  -8G+4=0  which  results  when 
in  (44)  with  k-|=0  is  equated  to  Qg  in  (22).  Thus,  G=2/3  is  the 

abscissa  of  the  point  in  the  QG-plane  in  Figure  8b  at  which  the 
flutter  boundary  and  the  divergence  (of  the  first  type)  boundary 
intersect. 
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It  is  obvious  from  (48)  that,  for  a double  pendulum  subjected 
to  a tangential  force,  sufficiently  large  damping  raises  the  value 
of  Qj  beyond  that  of  Qg^  as  given  in  (45).  However,  it  should  be 
remarked  that,  under  certain  circumstances,  even  very  large  damping 
may  be  of  a destabilizing  nature  when  a^l . For  example,  suppose 
that  <2=G=0,  so  that  the  difference  of  (49)  and  (45)  is 

AQ  = [11  + <1  - 3a(5+K^ )]/6(l -a).  (50) 

In  the  event  that  a=3/5  (for  the  sake  of  being  specific),  (50) 

becomes  AQ=(5-2k-|  )/6.  Hence,  for  sufficiently  large  Oi  Qd<.Qe2 

whenever  Ki<5/2  and  Q.  <Q.  whenever  k >5/2. 

I — o d I 


CONCLUSIONS 

By  means  of  example,  it  has  been  shown  that  a non-dissipative, 
quas i -dynami c system  composed  of  a certain  double  pendulum  subjected 
to  a circulatory  force  and  some  conservative  forces  can  lose 
stability  only  by  divergence.  Two  modes  of  onset  of  divergence  have 
been  identified:  one  associated  with  a vanishing  natural  frequency 

and  the  other  with  a natural  frequency  becoming  infinite.  Multiple 
regions  of  stability  and  instability  were  observed  in  apposite 
stability  maps.  However,  if  frictional  forces  are  acting  in  the 
hinges  of  the  double  pendulum,  the  secondary  stable  regions  disappear 
and  the  second  mode  of  divergence  is  replaced  by  loss  of  stability 
due  to  flutter.  Furthermore,  if  the  damping  forces  are  arbitrarily 
small,  the  value  of  the  critical  flutter  load  was  found  to  be 
always  less  than  the  critical  divergence  load  Qg^  for  the  corresponding 

undamped  system.  In  addition  to  the  results  reported  here,  still  other 

aspects  of  certain  quasi -dynamic  systems  have  been  discussed  in  refer- 
ences [7]-[10]. 
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Figure  1.  A double  pendulum. 
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Figure  2a.  Eigencurves  for 
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Figure  2b.  Eigencurves  for  o = 1/2. 
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Figure  2d.  Eigencurves  for  a = 3/2. 


27 


Figure  3.  Stability  map  for  a quasi -dynami c system  (k,=kp=G=0). 
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The  influence  of  weight  on 
a q ua s i -dynami c system  for 


the  stabi 1 i ty  of 
a = 0 (k]=k2=0). 
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Figure  4b.  The  influence  of  weight  on  the  stability  of 

a quasi-dynami c system  for  a = 3/4  (k]=k2=0). 
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Figure  4c.  The  influence  of  weight  on  the  stability  of 
a quas  i -dynami  c system  for  a = 1 (ki=l<2  = 0). 
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Figure  5a.  Stability  map  for  the  system  with  an 
end  constraint  for  a = 0 (G  = k-|=0). 
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Figure  5b.  Stability  map  for  the  system  with  an 
end  constraint  for  a = 1/4  (G  = k-|=0). 
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Figure  5c.  Stability  map  for  the  system  with  an 
end  constraint  for  a.  = 1 (G  = k]=0). 
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Figure  5d.  Stability  map  for  the  system  with  an 
end  constraint  for  a = 3/2  (G=ki=0). 
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Figure  6a.  Stability  map  for  the  system  with  a central 
constraint  for  a = 0 (G=k2=0). 
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Figure  6b.  Stability  map  for  the  system  with  a central 
constraint  for  a = 2/5  (G=k2=0). 
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Figure  6c.  Stability  map  for  the  system  with  a central 
constraint  for  a = 1/2  (G=k2=0). 
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Figure  6d.  Stability  map  for  the  system  with  a central 
constraint  for  a = 0.504  (6=l<2  = 0). 
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Figure  6e.  Stability  map  for  the  system  with 
constraint  for  a = 3/5  (G  = l<2  = 0). 
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Figure  7.  Stability  map  for  a slightly  damped  system  (k 
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Figure  8a, 


The  influence  of  weight  on  the  stability  of  a 
slightly  damped  system  for  a = 3/4  (ki=k2=0). 
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Figure  8b.  The  influence  of  weight  on  the  stability  of  a 
slightly  damped  system  for  a = 1 (k]=k2=0). 
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Figure  9a.  Stability  map  for  the  slightly  damped  system 
with  an  end  constraint  for  a = 1/4  (G=k]=0). 
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Figure  9b.  Stability  map  for  the  slightly  damped  system 
with  an  end  constraint  for  a = 1 (G  = k-|=0). 
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Figure  9c.  Stability  map  for  the  slightly  damped  system 
with  an  end  constraint  for  a = 3/2  (G=k]=0). 
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Figure  10a.  Stability  map  for  the  slightly  {lamped  system  with 
a central  constraint  for  a = 0.504  (G=k2=0). 
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Figure  10b.  Stability  map  for  the  slightly  damped  system  with 
a central  constraint  for  a =.3/5  (G=k2=0). 
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